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ABSTRACT 

We study the solutions of the equations of motion in the gauged (2+l)-dimensional 
nonlinear Schrodinger model. The contribution of Chern-Simons gauge fields leads to a 
significant decrease of the critical power of self-focusing. We also show that at appropriate 
boundary conditions in the considered model there exists a regime of turbulent motions 
of hydrodynamic type. 
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1. Introduction 



The nonlinear Schrddinger equation (NSE) is one of the basic models for nonlinear waves. 
The conventional field for NSE application is nonlinear optics Jl], |2[ wherein it describes 
propagation of wave beams in dispersive nonlinear media. The NSE appears also when 
one considers various nonlinear waves in hydrodynamics and plasma physics ||. A most 
important application area is the problem of detailed description of collapsing field distri- 
butions for the NSE with local cubic nonlinearity @, |J . In the case of the opposite sign 
of nonlinearity the NSE is used as the basic model f| of the low-dimension field theory 
to describe vortices in the problem of Bose condensation [[/]. 

Recent interest in problems with NSE solving in spatially 2D systems was associated 
with the account of the specific character of 2 + ID systems, which is manifested in 
equipping the NSE with a gauge field by replacing the usual derivatives by covariant ones. 
By that, the gauge field satisfies its own equation of motion (EM) with current, determined 
by the NSE solution. In the infrared limit, the main contribution to the EM of the gauge 
field in 2 + ID system is given by the Chern-Simons (CS) term within the system under 
consideration. The contribution of the gauge fields compensates in the Hamiltonian the 
contribution of nonlinearity at a certain relation of coupling constants. It leads to a soliton 
distribution of the fields found in Ref. ||. The nature of this phenomenon can easily be 
understood if one takes into account that CS term breaks T- and P-inversion symmetry 
of time and space. The chosen direction of the vector in the direction perpendicular to 
the plane can be considered as a chosen direction of rotation in the plane leading to the 
appearance of efficient repulsion. When this repulsion compensates the attraction, the 
Hamiltonian turns out to be limited at the bottom, and the CS solitons || correspond to 
its zero value. These field distributions are solutions of the self-dual equations ||. 

The results of the Refs. || [II], |TTJ stimulated a series of papers in this field. We would 
like to pay attention to some of them. In Ref. jnj the structures of field configurations 
were analyzed in full detail for a nonlinearity in NSE, which describes repulsion (in the 
absence of CS interaction) and takes into account the contribution of the non-zero mean 
value of particle number density. Considering the initial problem the authors of Ref. |13| 
concluded that at the most general initial condition the problems of EM for NSE with 
CS gauge fields correspond to the collapse regime. However, neither the spatial structure 
of the collapsing mode, nor the critical power for it (i.e. the number of particles in the 
mode) have been analyzed in this paper. The problem of exact integration of the model 
under consideration was analyzed in Ref. [fUR . The main result obtained is that the 
system cannot be integrated exactly excluding the following two cases: self-dual limit || , 
and the situation when the 2 + ID equations can be reduced to the 1 + ID equations. 
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An additional but still very important result of that paper was that the solitons of the 
gauged nonlinear Schrodinger equation (GNSE) have movable singularities on some curves 
in the two-dimensional plane. Detailed study of topological defects in low-dimensional 
systems have always been a key point for understanding dynamics of field distributions. 
The problem of so-called semi-local topological defects in the Chern-Simons-Higgs model 
was considered in a recent paper ]15| . 

The distribution of complex-valued functions $f(x,y,t) are defined on the manifold 
A4 which is multiconnected in spatial 2D systems. Therefore the fundamental homotopy 
group 7Ti(A4) determining analytical properties of function \l/ coincides with the braid 
group. There are actually several equivalent ways to reflect this fact in the theory. One of 
them is the Lagrangian approach that includes the effect of CS gauge fields into consid- 
eration. The CS term codes the existence and specific character of 2D point peculiarities 
contained in the Aharonov-Bohm gauge potentials within the long-wave description. One 
usually speaks of the long-range interaction represented by means of the CS gauge field as a 
statistical interaction between different field configurations. The different representations 
induce the different forms of field distribution. There is the so-called anyon representation 



T6| , [IT|| , when the gauge field in the explicit form is excluded from the Hamiltonian of 



the model thus providing the representation of " non- interacting" (by means of the gauge 
field) configurations of field ^(x, y,t). However, in this case the gauge field is proved 
to be included into the phase of function fy(x,y,t) that contains a cut in the complex 
plane, which provides multivaluedness of the function. The cut describes a string attached 
to the point defect. Thus, the so-called non-local topological defect separates sheets in 
the multi-sheet coverage of the 2D basic space. It is well known [lTj that the rep- 
resentation with explicit presence of gauge fields in the Lagrangian has a different form 
depending on the parity of the permutation group representation. The approach, based 
on the representation with fractional statistics of field configurations, is identical to the 
dynamic approach, when we are interested in the influence of a "statistical" gauge field. 

The existence of a gauge interaction has an exclusively topological character and is 
not associated with the quantum theory. This interaction, as a rule, was not taken into 
account when studying the classical dynamics of nonlinear models with a complex field 
in spatial 2D systems. Topological peculiarities certainly put additional restrictions to 
the quantization procedure in these systems []17|| . The role of the CS gauge interaction in 
this case is to take into account the vortex part of phase dynamics, which was not usually 
considered in classical systems when the 2 + ID NSE model was used. 

The purpose of this paper is the study of EM in the 2 + ID GNSE model. The main 
attention is given to the investigation of the structure of collapsing distribution of the 
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fields. Specifically we find, by means of numerical integration of the EM, the dependence 
of the critical power and efficient width of the zero-energy mode on coefficient k before 
the CS term. The limit k —>■ oo, when interaction with the gauge field is negligible, may 
be used as a test. In this case the known values of the power and the width are restored. 

If the phase of the field y, t) describes the longitudinal part in the gauge potential 
completely , evolution of field configurations is determined only by temporal dependence of 
the gauge field. We show that in this case the equations for the gauge field coincide with 
the EM of an ideal fluid. The effects of the manifestation of gauge field in classic systems 
with nontrivial topology, including the swimming motion at low Reynolds number within 
the 2 + ID hydrodynamics are well known [TJ| . A new feature is the fact that the basis 
for the 2D turbulence basing from the Euler equation in this case is chaotic dynamics of 
the CS gauge field. In this sense, the GNSE is a useful hydrodynamical tool |19| . 



One can see the following link between the dynamics of the CS fields and the problem 
of 2D turbulence. It is well known that CS action with appropriate boundary conditions 
is a way to classify conformal field theories PU| . The tools of the conformal field theory, 
in its turn, may be used |^T| to study the 2D turbulence. We show that within the 
model under consideration the connection between the dynamics of CS fields and 2D 
turbulence problem may be stated beyond the application of the conformal field theory. 
This dependence can be represented considering the evolution of closed current lines 
(loops) with the stochastization of lines near the points of link of the loops. The effect 
of contour links in terms of this paper reflects the effect of braiding world lines of the 
Aharonov-Bohm point singularities with formation of knots and links after projection 
of world lines onto the 2D space. Stochastization near the contour link points within 
formulation of 2D hydrodynamics of an ideal fluid in terms of contour variables |23|, |2^ 
was discovered in Ref. f24| . Such a stochastic behavior has a universal character. It is 



based on existence of the braid group and is closely connected with the arbitrary character 
of localization along the "time" axis of the point of world lines interlacing. Because of 
this the index k of the linking number proves to be a hidden parameter, which is not 
included into the Euler equations explicitly. 

The paper is organized as follows. The second section contains the EM for the CS 
gauge field and the GNSE for two different Ansatz corresponding to the goals of this paper. 
The third section is devoted to the numerical analysis of the problem. In conclusion, the 
results and open questions are discussed. 
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2. Equations of motion 

We consider a system with a Lagrangian density 

C = ^AadpA, + M*(d t + L4 )* i&m 2 + 1 1^| 4 • (1) 

The equations of motion (EM) have the form 

ify* = - iA) 2 ^ + A ^ - g\^\ 2 q , (2) 

[Vx A] ± = -i|*| 2 , (3) 

d t Ai + diA = ■ (4) 

Here g is the coupling constant and j = Jm\l/*(V — iA)^f is the current density. Hamil- 
tonian for Eq. (1) is 

H = X - j dh - tA)^\ 2 - g\^) , (5) 

where the potential A^ which is the auxiliary variable is expressed in terms of \^\ 2 in the 
following way 

A(r,t) = l|rfVG(r-r')|^| 2 (r',t) ! (6) 
A (r,t) = i|rfVG(r-r')j(r',t), (7) 



where the Green function G(r) 
satisfies the equation 



V x G(r) = -5 2 (r) (9) 

thus A^ is the solution of Eqs. (3) and (4). Since in Hamiltonian formulation the potentials 
are unambiguously presented by Eqs. (6) and (7), the gauge freedom 

Ap^A^-dtf, (10) 

tf-^tf (11) 
is fixed. This is achieved by the Coulomb gauge V • A = 0, with boundary conditions 

lim r 2 Ai(r, t) = J- £i -a-TV , (12) 

lim A (r,t) = 0. (13) 
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The choice of the boundary condition (12) is associated with the necessity to satisfy the 
integral representation of Gauss law (3) of CS dynamics 

$ = y>r[Vx A] ± = ~ J d 2 r\V\ 2 . (14) 

Here the magnetic flux $ and the number of particles 

N — ( d 2 r\m\ 2 (15) 



are conserved giving the global constrain $ = —\N. In the result of Eqs. (2)-(4), there 
exists naturally the continuity equation 

<9^| 2 + V-j = 0, (16) 

Let us use dimensionless fields and coordinates obtained by the following substitutions 

k 2 

q? = \hf 2 P e 1 *, A = -—w-d t tp, A x = -ku + d x tp, A y = -kv + d y tp, (17) 

' ' 2 t, v-itt- (is) 



fc|ife| \k\ \k 

The EM and the continuity equation, expressed by the new real functions p = p(x, y, t), 
u = u(x, y,t), v = v(x, y,t), w = w(x, y, t) have the form 

Pxx + Pyy = -2Cp 3 - pW + P [U 2 + V 2 ) , (19) 

u y -v x = -p 2 , (20) 

u t -w x = -2vp 2 , (21) 

vt-w y = 2up 2 , (22) 

p\ = 2 {{up 2 ) x + (vp 2 ) y ) (23) 

with the parameter C = g\k\ and notations u t = d t u etc. 
In the case of usual NSE 

id t m = -v 2 m - \m\ 2 m (24) 

substituting \I> = pe~ t<f ^ x ' y ^ we have 

Pxx + Pyy = -P 3 - PVt+P ({Vx) 2 + iVyf) , (25) 

P \ = 2 (( VxP 2 ) x + ^y P 2 )y) ■ (26) 

Comparing Eqs. (25), (26) with Eqs. (19) and (23), we pay attention to the following 
distinctions. First, due to gauge invariance, there are no derivatives of phase f in Eq. (19) 
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which exist in (25). Their role is played by gauge potentials. Therefore the evolution of 
the field p(x, y, t) is defined by the time derivatives of the functions u(x, y, t) and v(x, y, t) 
in Eqs. (20), (21). Unlike Eq. (25) the fields u and v are responsible for the transverse 
dynamics of the phase of the field \E r . Longitudinal dynamics of the phase is described by 
the zero component w(x, y, t) of the gauge potential which takes the place of the function 
ip t in Eq. (19). The function w(x,y,t) plays the role of a Lagrange multiplier permitting 
to take into account the restriction (20) of the Gauss law $ = — -rN locally. 

Second, the continuity equation (23) can be obtained excluding the function w from 
Eqs. (20) and (21) with the aid of Eq. (20). This remark is associated with the following 
problem. Let us assume in the Coulomb gauge V • A = — u x — v y + Aip = that the 
phase ip satisfies the equation Aip=0. Then the solution to the equation u x + v y = may 
be expressed by means of the function a(x, y, t) in the following way: 

u = dy, v = -a x . (27) 

In this case after replacing t by —It Eqs. (20) and (23) have the form 

a xx + a yy = -p 2 , (28) 

p 2 t + up 2 x + vp 2 y = 0. (29) 

The set of Eqs. (28) and (29) represents the "vorticity" form of Navier-Stokes equa- 
tions (Euler equations) for two-dimensional flows of ideal incompressible fluid where the 
function a(x,y,t) has the sense of the stream function. Note that hydrodynamic analo- 
gies have already been used for the solution of 1 + ID NSE problem p5|, |19fl . However, 
the exact proof that dynamics of CS gauge field in the frames of GNSE model (in the 
particular case of Coulomb gauge with Aip = 0) is equivalent to the two-dimensional EM 
of ideal incompressible fluid is given in this paper for the first time. The remarkable fact 
is that there is a close analogy between the states with the constant flux in the turbulence 
and CS anomaly [^TJ exposed by Eq. (28). 



It is useful to compare the gauge invariance of the model and the used Coulomb gauge 
at A(p = with canonical transformations and with area-preserving transformations. The 
infinitesimal area-preserving diffeomorphism which acts in the frames of CS theory has 
the form 

£i-*€i+Xi> d iXi = o , (30) 

where £j = (x, y), x% — (^i> ^2)- The general solution of the equation di\i = is the sum 
of two terms 

bi 

Xi = £ijdja(£) + J2 °kXi , (31) 



k=l 
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where the second term describes the finite number (given by the first Betti number b\) 
of harmonic forms on the two-dimensional phase space (A x , A y ) of CS theory. Diffeo- 
morphisms which resulted from the first term in Eq. (31) are nothing but canonical 
transformations p6| . In the case of phase space which is actually a torus there are also 



two global translations p& = —idk (for torus the first Betti number b\ — 2). Just in this 
case the phase ip(x, y, t), satisfying the equation Aip = 0, is the linear function ip = ax+by. 
From the viewpoint of NSE this corresponds to the constant direction of ray propagation 
assigned by the vector n ~ (a, b). In the general case of the phase space with arbitrary 
topology it is invalid and the phase (p(x,y,t) not satisfying the equation Aip = gives 
rise to the "additional" longitudinal contribution to the potentials u(x, y, t) and v (x, y, t). 

Let us consider for example the case when the Ansatz for the field ^/(x,y,t) corresponds 
to the generalized lens transformation || |13| 

vl/ (r, t) = ?^ exp (-ib{r)e/2 + l \r) . (32) 

Here £ = r/g(r), r = / du [f(u)\ , 6(r) = —f t f = —g T 9- The gauge potentials transform 

o 

M, at such substitution, follows 

A(r,t)-[g(r)]- 1 A(Cr), (33) 

A (r, t) - [g(r)]- 2 [MC r) - 6(r)CA(C r)] (34) 

the relations (6), (7), being preserved, where the function p = |$|. After these transfor- 
mations Eq. (2) changes its form 

id T § + (PC ~ A)$ = --(V - 2A) 2 $ + - #|$| 2 $ , (35) 

because the function /3(r) = {b 2 + b T )/2 = —f 3 f tt /2 in the case of <p(x,y,t) ~ b(x 2 + y 2 ) 
and b(t) ^ t — t does not equal zero. If we are interested in collapsing solutions with 
[p7| , |i~3H }' 2 {t) ~ (to — t)/ln [In (to — t)], the structure of the self-similar nonlinear core |13| 
of the solution is described by the solutions of the following equation 

-A$ = --(V-zA) 2 $ + A $-#|$| 2 $. (36) 

In the next section, by numerical calculation, we find the zero-energy localized ground 
state of GNSE (36). We show the dependence of its effective width on the parameter 
C — g\k\ as well as the form of the functions u, v, w. 
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3. Solution structure 



For the numerical analysis of the solutions of Eq. (36) we use the method of the stabilizing 



multiplier [28]. The iteration approach for Eq. (36), which differs from Eq. (19) by the 
additional term — A$ in the LHS has the form 

= M n F- x (G(p)F (-2C$l + j<S> n (u 2 + v 2 - w) n )) , (37) 

M = ( fdMF* n y V m) 

n \fdZpG(p)F<S> n F(-2C& n + j$ n (u* + vZ-w) n )) • 1 ] 

Here F (F^ 1 ) are the operators of direct (inverse) Fourier transform, G(p) = — {p 2 + A) . 
The multiplier j — 1 or j — if we take into account the nonlinear contribution of 
gauge fields in Eq.(36) or neglect it. In the case j = the usual normalization in NSE 
corresponds to C = 1/2. Without the restriction of generality we shall suppose below 
that A = 1. 

We should choose the exponent a in the stabilizing multiplier M n comparing the 
degrees of homogeneity of terms in LHS and RHS of Eq. (36) proceeding from the 
requirements that M n — > 1 at n —>■ oo. Without the term $(u 2 + v 2 — w) the exponent 
a equals 3/2. New features of our problem are that the nonlinearity in Eq. (36) has the 
polynomial character of the type — 2C$ 3 +6$ 5 because both of the terms and §(u 2 +v 2 ) 
are proportional to $ 5 . Therefore, for the convergence of the iteration approach a should 
belong to the range 5/4 < a < 3/2. In the simulation of the present paper we have used 
the value a = 3/2 which gives rapidly the value M n = 1 of the stabilizing multiplier. 
We have used the distribution of the form $(Cx, C y ) = (VT?/ 71 ") ex P( — 7Ca; 2 — d(y 2 ) as the 
initial field configurations with arbitrary constants 7 and 5. In our calculations we have 
obtained rapidly the isotropic solutions. 

To regularize the integrals (6) and (7) which diverge logarithmically in coincident 

points at numerical calculations of gauge potentials u and v we substituted r 2 — >• r 2 + e 2 

in the expression for Green function (8). In the momentum space this corresponds to 

00 

the substitution of d 3 pf(ep) with f(ep) = J dmexp(~y/m 2 + e 2 p 2 ) for d 3 p. The factor 



f(ep) at ep 3> 1 decreases exponentially cutting off all momentum integrals. However, 
the infrared region remains the same, because at ep 1 f(ep) = 1. In our calculations 
the cut-off radius e which has the sense of the thickness of vortex core, was equal to 10~ 2 . 

Simulation was performed on the square lattice with linear sizes L x = L y = 12. The 
maximum number of lattice sites was limited by the value n = n x n y = 128 • 128. To test 
our approach we used the solution of EM (36) with = (j = 0) and with C = 1/2 
which gives the well-known value N = 11.703, as well as the solution of a self-dual equation 
A In p = —p 2 H , when w = —p 2 , u — d y In p, v — —d x In p. 
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The Figures 1-3 show the configurations of the fields p, u and w for the specific value 
of the parameter C = 4. We may obtain the form of the function f(Cx,Cy) using the 
relation v(( x , Q) = -u(( y ,Cx)- 

Using the obtained function p we computed the dependencies of the critical power N 
(the particle number) and the effective width < R 2 >= iV -1 / <i 2 CC 2 p(0 on the parameter 
C. The results of calculations are given in the Table 1. 



Table 1. 



J 


c 


N 


< R 2 > 





0.5 


11.703 


1.2607 


1 


3 


2.9216 


1.2464 


1 


5 


1.2825 


1.2579 


1 


10 


0.5973 


1.2600 


1 


100 


5.8528-10" 2 


1.26066 


1 


1000 


5.8516-10" 3 


1.26066 



4. Conclusion 

In this paper we have studied systematically the influence of the CS gauge field, reflecting 
the specific feature of the dimensionality of our problem, on field configurations in the 
GNSE model. Here we summarize some new results. 

It is seen from Eqs. (15)-(19) that if we neglect CS gauge fields (j = in Eq. (37)) the 
dependence of the particle number iV on the parameter C = g\k\ has the form N = Nq/C. 
This dependence is shown in Fig. 4 by the dotted line. It follows from the results shown in 
the first line of Table 1 that Nq = 5.585. The contribution of CS gauge fields (j=l in the 
Table 1) leads to the sharp decrease of values of N. In particular iVj =1 (3)/iVj =0 (0.5) ~ 
0.25. The effective width < R 2 > changes slightly. The calculated dependence N(C) at 
C > 3 is given in Fig.4. 

As it should be expected at a fixed value of the parameter C in the region C > 3, 
Nj=i(C) is always greater than Nj =Q (C), because the CS gauge fields describe effective 
repulsion. In the range 1 < C < 3 we could not perform calculations in the framework 
of the method used, due to the breaking of convergence of the iteration method (37), 
(38). We cannot explain this phenomenon using the theory presented in this paper. Note 
however, that the values C = 1 and C = 2 at g = 1 correspond to the discrete values 
k = 1 and k = 2 which describe the field configurations in the self-dual limit k = 1 || 
and the one-half fermion limit k = 2, which qualitatively differ from our case. The classic 
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limit of the considered theory corresponds to the case k — > oo, when the gauge field splits 
off from the field V(x, y, t). For Eq. (19) and for the value N(C) the limit C > 1 denotes 
that N(C) -> 0. 

Strong Langmuir turbulence in plasmas is usually described by the solutions of the 
NSE (Eq. 24). It is assumed that a cascade of the randomly distributed self-similar 
collapsing fields is generated. In this paper we show that the specific features of spatial 
two-dimensional systems may lead to the traditional picture of turbulence associated with 
Euler equations. However, for the hydrodynamical mechanism of turbulence (HMT) to 
be involved, it is necessary the existence of a linear profile of the phase <p(x, y) (= ax + by) 
in each mode. This implies that the nonlinear contributions (by x and y) in the temporal 
evolution of the phase are small. 

One of the media in which the HMT can play a role is an optical medium with random 
inhomogeneous guiding surfaces. Reflecting from the surfaces, wave fronts acquire random 
directions of propagation. For media with weak Kerr nonlinearity, the nonlinear phase 
disturbance from adjacent points will not be important. Application of the HMT model 
suggested above requires separate consideration and will be presented elsewhere. 

In conclusion, we studied numerically the structure of the collapsing mode in GNSE 
model, observed the effect of strong reduction of the critical power N in spatial two- 
dimensional systems as compared to the traditional values, and showed that in the case 
of appropriate boundary conditions the phenomenon of collapse inhibits the development 
of turbulence according to the hydrodynamic scenario. 
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Figure captions 

FIG. 1. Plot of the function p(( x , 0) = p(0, Q. 
FIG. 2. Plot of the function u(( x ,0). 
FIG. 3. Plot of w(( x , 0) = w(0, Q. 

Fig. 4. Number of particles N as a function of the parameter C without 
taking into account the gauge field (dotted line) and with the gauge 
field (solid line). The point denotes the value iV(0.5) = 11.703. 
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